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ABSTRACT 

cRNA microarray and real-time PCR (qPCR) studies identified a number of Estrogen Signaling Pathway 
(ESP) genes (GTF2H2, NCOR1, TAF9, NRAS, NRIP1, POLR2A, DDX5, NCOA3) and Cell Cycle Pathway 
genes (CCNA2, PCNA, CDKN1A, CCND1, PLK1) in MCF-7 breast cancer cells that are regulated by the 
bioflavonoid luteolin. Chromatin immunoprecipitation (ChIP) studies revealed that luteolin modified histone 
H4 acetylation at the PLK-1 promoter suggesting that this bioflavonoid controls gene transcription via an 
epigenetic mechanism involving histone H4 acetylation. These findings are consistent with the anti-estrogenic 
and anti-proliferative properties of luteolin in normal and malignant cells. {IntJBiomed Sci 2011; 7 (2): 101-111) 
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INTRODUCTION 

The bioflavonoid luteolin demonstrates anti-estrogenic 
and/or anti-proliferative properties in a variety of model 
systems including the rodent uterus and prostate and hu- 
man breast and prostate cancer (1-3). The regulatory ef- 
fects of luteolin on normal and malignant cell prolifera- 
tion likely involves its high affinity binding interactions 
with nuclear type II [ 3 H]estradiol binding sites that were 
recently identified as histone H4 (4-6). That type II sites 
represent ligand binding to histone H4 is consistent with 
the ubiquitous nature of type II binding sites in mam- 
malian cells and the abilities of type II site ligands with 
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anti-estrogenic and anti-proliferative activities to inhibit 
the proliferation breast, prostate, colorectal, leukemia, me- 
ningioma, and ovarian cancer cells (1,2, 7-15). The abili- 
ties of luteolin and related ligands to regulate normal and 
malignant cell proliferation are attributed to their capaci- 
ties to mimic methyl-p-hydroxyphenyllactate (MeHPLA) 
as cell growth regulatory agents. MeHPLA was identified 
as the endogenous ligand for nuclear type II sites (16), and 
an esterase-induced deficiency of MeHPLA in malignant 
tissues and cells directly correlates with the loss of regula- 
tory control (16, 17). Consequently, esterase stable MeH- 
PLA-related compounds such as luteolin are capable of 
blocking estrogenic response in the uterus and the growth 
of the prostate and breast and prostate cancer cell prolif- 
eration in vitro and in vivo (1, 2, 7, 8, 17, 18). 

Nearly three decades of work in our lab indicate that 
MeHPLA and nuclear type II sites are components of a 
classical cell growth regulatory pathway in mamma- 
lian cells. That nuclear type II sites are a ligand-binding 
domain on histone H4 further suggests that ligands for 
this binding site may control specific gene transcription 
through an epigenetic mechanism influencing histone 
acetylation and/or deacetylation events. Luteolin appears 
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to regulate the transcription of a number of genes in the 
epidermal growth factor signaling pathway (EGFSP) and 
cell cycle pathway (CCP) in human prostate cancer cells 
through such an epigenetic mechanism involving histone 
H4 acetylation (3, 19). The microarray, qPCR and ChIP 
studies in the present manuscript further suggest that the 
antiestrogenic and anti-proliferative properties of luteolin 
in MCF-7 breast cancer cells may be mediated through 
the regulation of genes in the estrogen signaling pathway 
(ESP) and CCP's. ChIP studies suggest that this regulation 
involves histone H4 acetylation. 

MATERIALS AND METHODS 

MCF-7 Cells and Growth Conditions 

Stock cultures of MCF-7 cells were grown in T-150 
flasks in DMEM containing 5% FBS (20). Luteolin inhib- 
its the proliferation of our ER(+) MCF-7 cells (Figure 1) 
which contain ERa (but not ERP) even though the flavo- 
noid does not bind to either form of the ER (21). For the 
microarray studies, MCF-7 cells were grown in phenol- 
red free DMEM supplemented with 5% charcoal-stripped 
fetal calf serum (CSFCS) for 48 hours. The cells were 
then re-fed with fresh phenol-red free DMEM. At this 
time, triplicate T-75 flasks for each treatment group were 
treated for 5 hours with estradiol- 17 (3 (5 nM) or luteolin 
(17.5 uM) alone or in combination added to the medium 
in 2 uL of ethanol (vehicle controls). This level of luteolin 
inhibits MCF-7 cell proliferation by about 50% at 48 hours 
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Figure 1. Luteolin Inhibition of Cell Proliferation. MCF-7 
cells were treated for 48 hours with 2 )iL of ethanol or the indi- 
cated concentrations of luteolin in 2 uL of ethanol and attached 
cells (mean ± SEM for 12 wells) were determined by cell counts 
(Trypan Blue exclusion) and expressed as a % of the ethanol 
controls. 



following treatment (Figure 1) and appears cytostatic but 
non toxic (16, 21). Cell number was determined by Trypan 
Blue dye exclusion as previously described (3). For valida- 
tion of the microarray data, RNA was prepared for qPCR 
analysis from MCF- cells grown as described. Triplicate 
10 cm dishes of MCF-7 cells were treated with estradiol 
(5 nM) or luteolin (17.5 uM) alone or in combination for 5 
hours prior to harvesting and RNA preparation (19). 

RNA Preparation for Microarray or qPCR Studies 

MCF-7 cells from the various treatment groups of 
flasks or dishes were washed with PBS and collected with 
0.25% trypsin-EDTA (19, 20). Approximately 5.0 x 10 6 
cells from each flask or plate were collected by centrifu- 
gation (2000 rpm x 5 min) in RNAse/DNAse free tubes 
and resusp ended in 1 mL of PBS plus 4 mL of RNAlater 
(Qiagen) and stored at -20°C. The frozen cells were thawed 
on ice, collected by centrifugation and lysed in 0.6 mL of 
RTL buffer (Qiagen) containing P-mercaptoethanol. The 
lysed cells were passed Qiashredders (Qiagen) and the 
passthrough diluted with an equal volume of 70% EtOH 
and loaded onto RNeasy spin columns (Qiagen). The col- 
umns were washed with RW1 followed by RNAse-free 
DNAse digestion and the purified total RNA was eluted 
with 50 uL of RNAse free water as per the manufacturer's 
instructions. RNA integrity was confirmed on an Agilent 
2100 Bioanalyzer in the Baylor College of Medicine Mi- 
croarray Core Laboratory (21). 

Microarray Analyses 

RNA samples were subjected to oligo(deoxythymidine)- 
Reverse Transcription, in vitro transcription and biotin- 
labeling of cRNA (Enzo Biochem, Farmingdale, NY). 
The fragmented, biotin-labeled cRNA was hybridized 
to Human Genome U133 plus 2.0 oligonucleotide arrays 
(Affymetrix, Santa Clara, CA) containing approximately 
54,000 probe sets (38,500 genes) as described (19, 20, 22). 
Each transcript is represented on the chip as 1 1 probe pairs, 
each pair containing a perfect match and mismatch. Micro- 
arrays were stained with streptavidin antibody and strepta- 
vidin-phycoerythrin in an Affymetrix Fluidics station and 
scanned at 3 mm with a GeneArray scanner (Affymetrix). 
Experiments were performed in triplicate with indepen- 
dent pools of cRNA from EtOH controls, estradiol, luteo- 
lin or estradiol plus luteolin treated MCF-7 cells using 12 
separate microarray chips (three chips for each group). Fol- 
lowing low-level quantification of the scanned data using 
GeneChip Operating System (GCOS, Affymetrix), data 
were further analyzed with dChip 2006 (Harvard) to adjust 
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the arrays to a common baseline and to estimate expres- 
sion using the PM-only model (23, 24). The 12 Genechips 
were normalized to the same baseline (Genechip with the 
median average intensity; in this case ethanol control Chip 
#3) and modeled together (25). Data quality was reviewed 
using preset call rates from GCOS (average 33.63%: range 
28.9% to 41.7%), ratios of 3' to 5' glyceraldehyde 3-phos- 
phate dehydrogenase probe sets from GCOS (average 1.29: 
range 0.96 to 1.85) and array outlier rates from dChip (aver- 
age 0.33%: range 0.0% to 0.922%). Differentially expressed 
genes in the estradiol plus luteolin treated group relative to 
those in the estradiol treated group, were selected using a 
two-sample comparison with a lower boundary 90% confi- 
dence interval, fold change greater than 1.2 and a value dif- 
ference between group means of >50. The medium number 
of detected genes in 50 permuted samples was used to es- 
timate the false discovery rate. Similar comparisons were 
made for the EtOH controls vs. estradiol treatment group to 
identify the estradiol-regulated genes. The luteolin-regu- 
lated genes identified by comparison of the expression pat- 
terns in the estradiol vs. estradiol plus luteolin treated cells 
were essentially indistinguishable from those identified 
by comparing the expression in luteolin treated cells rela- 
tive to EtOH controls. Utilization of these four treatment 
groups generated valid comparisons of the gene expres- 
sion profiles for all groups. The microarray data were ana- 
lyzed with Ingenuity Pathway Analysis Software (Ingenu- 
ity Software, Redwood City, CA) or GenMAPP Software 
(Gene Map Annotator and Pathway Profiler, Version 2.1, 
Gladstone Institutes, University of CA at San Francisco) to 
identify differentially expressed gene pathways regulated 
by estradiol or luteolin according to Gene Ontology func- 
tion (3, 19, 20). 

Real-time Quantitative Polymerase Chain Reaction 
(qPCR) 

Luteolin effects on gene expression identified by the mi- 
croarray studies were subjected to qPCR analysis for con- 
firmation. cDNA for qPCR was generated from the RNA 
samples with the iScript cDNA synthesis Kit from BioRad. 
Commercially available, pre-validated primers (Qiagen) for 
the ESP and CCP genes were used in MyiQ SYBR Green 
Supermix and quantified on a MyiQ Single Color Real- 
Time PCR Detection System using MyiQ Optical System 
Software, Version 2.0 (Bio-Rad). Validation of each primer 
pair for the ESP and CCP genes was by generating stan- 
dard serial dilution and melt curves on cDNA from MCF- 
7 cells. Reaction efficiencies of 90-110% and correlation 
coefficients of >0.995 were obtained. Melt curves demon- 



strated a single reaction product with an appropriate melt- 
ing temperature confirming that primer dimerization was 
not contributing significantly to the signal. Results from 
quadruplicate qPCR runs on triplicate pools of RNA were 
normalized to 18 S RNA for each of the various treatment 
groups. Products of the optimized reactions were also an- 
alyzed by agarose gel electrophoresis to confirm that the 
amplicon size corresponded to the data provided by Qia- 
gen for each primer pair. Normalized qPCR data for each 
gene studied from estradiol plus luteolin treated cells were 
expressed as a percent of that observed in estradiol treated 
cells (controls for comparison to estradiol plus luteolin 
group) and analyzed by ANOVA and a suitable test on the 
treatment means utilizing Instat (Graphpad Software, Inc.) 
as previously described (3, 19, 20). These studies defined 
luteolin effects on the expression of each of the selected 
ESP or CCP genes in estradiol plus luteolin treated MCF-7 
cells relative to the estradiol-treated controls. 

Chromatin Immunoprecipitation of PLK-1 Gene 
Promoter 

Following the treatment of MCF-7 cells with estradiol 
or luteolin (alone or in combination) as described above, the 
monolayers were fixed with formaldehyde to cross-link pro- 
tein to DNA (19, 20). The chromatin was then digested (10- 
25 minutes) with Enzymatic Shearing Cocktail (ChlP-IT 
Kit, Active Motif) to generate DNA fragments with an av- 
erage length of about 200 base pairs. Aliquots of the digest- 
ed chromatin were pre-cleared with protein-A-sepharose 
and incubated with (+Ab) or without (-Ab) anti-acetylated 
(Lys 5, 8, 12, 16) histone H4 antibody. Following incubation 
with the protein-A-sepharose beads and cross-link rever- 
sal, the immunoprecipitated DNA was de-proteinized and 
analyzed by PCR and qPCR. Forward and reverse primers 
for the PLK-1 promoter (-262 to -80 bp) were 5'-GGTTTG- 
GTTTCCCAGGCTAT-3 ' and 5'-GCTGGGAACGTTA- 
CAAAAGG-3', respectively as described (26). 

RESULTS 

D-Chip Analysis of the Microarray Data 

Luteolin inhibition of MCF-7 cell proliferation is in- 
dependent of ERa or ER(3 binding interactions and this 
cell line is a well-accepted model for studying the inhibi- 
tory effects of type II site ligands on gene expression and 
cell proliferation (21). On the basis of the D-Chip analysis, 
genes regulated by luteolin in estradiol-stimulated MCF-7 
cells are shown in Table 1. Luteolin treatment of estradiol 
stimulated MCF-7 cells (estradiol + luteolin group versus 
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estradiol alone) down regulated 840 genes and up regu- 
lated 772 genes. The sum of regulated genes (1612) rep- 
resented 2.9% of 54000 probe sets (38,500 genes) on the 
Human Genome Affymetrix U-133 plus 2.0 chip. These 
genes were selected on the basis of a 90% confidence in- 
terval and a minimal fold change of 1.2. Thus, luteolin was 
highly selective in the regulation of specific genes in MCF- 
7 breast cancer cells 

Ingenuity Analysis of Microarray Data 

To identify biochemical pathways regulated by luteo- 
lin, the microarray data were first subjected to Ingenuity 
Pathway Analysis (Ingenuity Software, Redwood City, 



Table 1. Microarray Analysis of Luteolin Effects on 
Gene Expression in MCF-7 Cells 



Down 


% 


Up 


% 


Total 


% 


840 


1.54 


772 


1.41 


1612 


2.95 



MCF-7 cells were treated for 5 hours with estradiol (5 nM) or estra- 
diol (5 nM) plus luteolin (17.5 |j,M). Cells were harvested for RNA 
preparation and microarray analyses. Luteolin regulated genes 
were determinate by dChip comparison of the expression patterns 
between the cells treated with estradiol alone or estradiol plus lu- 
teolin. 



CA). In estradiol stimulated MCF-7 cells, luteolin sig- 
nificantly modulated the expression of genes in 15 out of 
112 canonical pathways identified with Ingenuity. Genes 
involved in purine metabolism (not shown) and the ESP 
(Table 2 and Figure 2) displayed the most pronounced re- 
sponse. Please note that all of the genes listed in Tables 2 
and 3 were components of the ESP or CCP, however, not 
all of these genes were depicted in the Ingenuity or Gen- 
MAPP pathway diagrams (Figures 2 and 4) which were 
down-sized due to space limitations. Key genes focused 
on in this manuscript are depicted in Figures 2 or 4. 

Although type II site ligands like luteolin do not bind 
to ERa or ER(3 (21), they are anti-estrogenic in the rat 
uterus (1, 2, 7, 18), and might be expected to modulate 
genes in the ESP. The data (Figure 2 and Table 2) show 
that within the ESP, estradiol stimulated the expression of 
4 genes including ESR1 (estrogen receptor 1; ERa), GT- 
F2H2 (general transcription factor II H, polypeptide 2, 44 
kDa), NCOR2 (nuclear receptor co-repressor 2) and PGR 
(progesterone receptor) and inhibited the expression of 6 
genes including CREBBP (CREB binding protein; Rubin- 
stein-Taybi syndrome), P300 (EiA binding protein p300), 
NCOA3 (nuclear receptor coactivator 3), POLR2J2 (DNA 
directed RNA polymerase II polypeptide J-related gene), 



Figure 2. Luteolin Effects 
on Estrogen Signaling 
Pathway (ESP) Genes. 

Ingenuity Analysis ofmicro- 
array data from MCF-7 cells 
treated with estradiol (con- 
trols) or estradiol + luteolin 
for 5 hours. Genes in this 
pathway were significantly 
(p<0.01) upregulated by 
luteolin (light blue sym- 
bols) or estradiol (shaded 
symbols) or significantly 
(p<0.01) downregulated by 
luteolin (gray symbols) or 
estradiol (dotted symbols). 
These data are numerically 
summarized in Table 2. 
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Table 2. Genes in the ESP Whose Expression was Significantly Changed by Estradiol (5 nM) or 
luteolin (17.5 uM) Treatment of MCF-7 Cells for 5 Hours. 


Gene 


Estradiol Regulated Genes 


Fold Change 


Entrez Gene 


CREBBP 


CREB binding protein (Rubinstein-Taybi syndrome) 


-1.57 


1387 


EP300 


E1A binding protein p300 


-1.89 


2033 


ESR1 


estrogen receptor 1 


1.75 


2099 


GTF2H2 a 


general transcription factor IIH, polypeptide 2, 44kDa 


1.68 


2966 


NCOA3 a 


nuclear receptor coactivator 3 


-1.65 


8202 


NCOR2 


nuclear receptor co-repressor 2 


1.76 


9612 


PGR 


progesterone receptor 


2.34 


5241 


POLR2J2 


DNA directed RNA polymerase II polypeptide J-related gene 


-1.69 


246721 


SMARCA4 


SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily 
a, member 4 


-1.45 


6597 


TRRAP 


transformation/transcription domain-associated protein 


-1.68 


8295 


Gene 


Luteolin Regulated Genes 


Fold Change 


Entrez Gene 


CTBP2 


C-terminal binding protein 2 


-1.43 


1488 


DDX5 


DEAD (Asp-Glu-Ala-Asp) box polypeptide 5 


-1.29 


1655 


GRB2 


growth factor receptor-bound protein 2 


-1.55 


2885 


GTF2H2 a 


general transcription factor IIH, polypeptide 2, 44kDa 


-1.43 


2966 


HNRPD 


heterogeneous nuclear ribonucleoprotein D (AU-rich element RNA binding protein 1, 
37kDa) 


-2.26 


3184 


NCOA3 a 


nuclear receptor coactivator 3 


-1.84 


8202 


NCOR1 


nuclear receptor co-repressor 1 


-1.52 


9611 


NRAS 


neuroblastoma RAS viral (v-ras) oncogene homolog 


-1.36 


4893 


NRIP1 


nuclear receptor interacting protein 1 


-1.35 


8204 


PCK2 


phosphoenolpyruvate carboxykinase 2 (mitochondrial) 


-1.65 


5106 


POLR2A 


polymerase (RNA) II (DNA directed) polypeptide A, 220kDa 


-1.58 


5430 


TAF5L 


TAF5-like RNA polymerase II, p300/CBP-associated factor (PCAF)-associated factor, 
65kDa 


1.48 


27097 


TAF6L 


TAF6-like RNA polymerase II, p300/CBP-associated factor (PCAF)-associated factor, 
65kDa 


1.9 


10629 


TAF9 


TAF9 RNA polymerase II, TATA box binding protein (TBP)-associated factor, 32kDa 


1.42 


6880 



a Modulated by both Estradiol and Luteolin. 



SMARCA4 (SWI/SNF related, matrix associated, actin 
dependent regulator of chromatin, subfamily a, member 
4) and TRRAP (transformation/transcription domain-as- 
sociated protein). 

In comparison, luteolin stimulated the expression of 
3 genes including TAF5L (TAF5-like RNA polymerase 
II, p300/CBP-associated factor (PCAF)-associated factor, 
65kDa), TAF6L (TAF6-line RNA polymerase II, p300/ 
CBP-associated factor (PCAF)-associated factor, 65kDa) 
and TAF9 (TAF9 RNA polymerase II, TATA box binding 
protein (TBP)-associated factor, 32kDa) and inhibited the 
expression of 11 genes including CTBP2 (C-terminal bind- 



ing protein 2), DDX5 (DEAD (Asp-Glu-Ala-Asp) box poly- 
peptide 5), GRB2 (growth factor receptor-bound protein 2), 
GTF2H2 (general transcription factor IIH, polypeptide 2, 
44kDa), FFNRP (heterogeneous nuclear ribonucleoprotein D 
(AU-rich element RNA binding protein 1, 37kDa), NCOA3 
(nuclear receptor coactivator 3), NCOR1 (nuclear recep- 
tor co-repressor 1), RAS (neuroblastoma RAS viral (v-ras) 
oncogene homolog), NRIP1 (nuclear receptor interacting 
protein 1), PCK2 (phosphoenolpyruvate carboxykinase 2 
(mitochondria)) and POLR2A (polymerase (RNA) II (DNA 
directed) polypeptide A, 220kDa). Only two genes (GTF2H2 
and NCOA3) were affected by both estradiol and luteolin. 
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Thus, the response of ESP genes to luteolin is specific and 
separate from genes regulated by estradiol in these cells. 

QPCR Analysis of ESP Genes 

To validate the microarray data (Figure 2, Table 2), trip- 
licate pools of RNA from MCF-7 cells treated with estradiol 
or estradiol + luteolin for 5 hours were converted to cDNA 
using the iScript cDNA Synthesis Kit (Bio-Rad, Hercules, 
CA) and the expression of L-regulated genes was confirmed 
by qPCR (19-21). Pre-validated primers for TAF9, DDX5, 
GRB2, GTF2H2, NCOA3, NCOR1, RAS, NRIP1, POL- 
R2A (purchased from Qiagen) were used for qPCR in MyiQ 
SYBR Green Supermix (Bio-Rad) by routine procedures in 
our lab (3, 19, 20). Results from qPCR on RNA from 3 sepa- 
rate experiments consisting of estradiol versus estradiol + 
luteolin treated cells were normalized to 18S RNA and ana- 
lyzed statistically by Instat (Graphpad Software, Inc.) with 
a suitable test on the treatment means (3, 19, 20). The qPCR 
data (Figure 3) show that luteolin stimulated TAF9 gene 
expression, as indicated by the microarray analysis (Table 



Table 3. Genes in the CCP Whose Expression was Significantly Changed by Estradiol (5 nM) or 
Luteolin (17.5 uM) Treatment of MCF-7 Cells for 5 Hr 


Gene 


Estradiol Regulated Genes 


Fold Change 


Entrez Gene 


CDC2 a 


Cell division cycle 2, Gl to S and G2 to M 


1.38 


983 


CDC20" 


CDC20 cell division cycle 20 homolog (S. cerevisiae) 


-1.45 


991 


CDC25B 


Cell division cycle 25B 


-1.21 


994 


CDKN1B 


Cyclin-dependent kinase inhibitor IB (p27, Kipl) 


-1.39 


1027 


GSK3b a 


Glycogen synthase kinase 3 beta 


-1.44 


2932 


MAD2L1 


MAD2 mitotic arrest deficient-like 1 (yeast) 


1.75 


4085 


PCNA a 


Proliferating cell nuclear antigen 


1.23 


5111 


PRKDC 


Protein kinase, DNA-activated, catalytic polypeptide 


-1.64 


5591 


Gene 


Luteolin Regulated Genes 






ASK 


Activator of S phase kinase 


1.43 


10926 


CCNA2 


Cyclin A2 


-1.45 


890 


CDC2 a 


Cell division cycle 2, Gl to S and G2 to M 


-1.27 


983 


CDC20" 


CDC20 cell division cycle 20 homolog 


-1.5 


991 


CKDN1A 


Cyclin-dependent kinase inhibitor 1A (p21) 


1.67 


1026 


GSK3b a 


Glycogen synthase kinase 3 beta 


-2.27 


2932 


PCNA a 


Proliferating cell nuclear antigen 


-1.35 


5111 


PLK1 


Polo-Like kinase 1 


-1.66 


5347 


PTTG1 


Pituitary tumor-transforming 1 


-1.29 


9232 



'Modulated by both Estradiol and Luteolin. 



2) and inhibited the expression of GTF2H2, NCOR1, RAS, 
NRIP1, POLR2A, DDX5 and NCOA3. 




GRB2 GTF2H2 NCOR1 TAF9 NRAS NRIP1 POLR2A DDX5 NCOA3 

Genes 

Figure 3. Luteolin Effects on ESP Genes. MCF-7 cells were 
treated for 5 hours with estradiol (5 nM) or estradiol + luteolin 
(17.5 uM) in 2 uL EtOH. RNA was isolated and gene expres- 
sion assessed by real-time PCR (qPCR). The results represent 
the mean ± SEM for three independent RNA sets normalized to 
18S RNA. *p<0.0001. 
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GenM APP Analysis of Microarray Data 

Analysis of the D-Chip data with GenMAPP also dem- 
onstrated the genes in the CCP were significantly (/?<0.001) 
affected by luteolin treatment. Of the 9 genes in the cell 
cycle pathway affected by luteolin (Table 3 and Figure 4), 
the expression of only CDKN1A (cyclin-dependent kinase 
inhibitor 1A or p21) ASK (activator of S phase kinase) was 
increased. The remaining 7 genes in this pathway were in- 
hibited. Interestingly, estradiol and luteolin both inhibited 
the expression of CDC20 (CDC20 cell division cycle 20 
homolog) and GSK3b (glycogen synthase kinase 3 beta) 
and these two compounds had opposite effects on CDC2 
(cell division cycle 2, Gl to S and G2 to M) and PCNA 
(proliferating cell nuclear antigen). 

QPCR Analysis of CCP Genes 

To validate the microarray data identifying the CCP 
genes as luteolin targets, a few genes in this pathway 
were selected for further analyses by qPCR (Figure 5). 
These included CDKN1A, cyclin A2 (CCNA2), PLK1 
and proliferating cell nuclear antigen (PCNA). Luteo- 
lin stimulated CDKN1A and inhibited the expression of 
CCNA2, PCNA and PLK1. Cyclin Dl (CCND1) was also 
examined because this gene is regulated by luteolin in es- 



tradiol- stimulated MCF-7 cells (21). For some unknown 
reason, CCND1 was not included in the GenMAPP Dia- 
gram (Figure 4), even though this is a well- known cell 
cycle gene (21, 27) whose expression was inhibited (Fig- 
ure 5) by luteolin. 

Chromatin Immunoprecipitation Studies 

Since luteolin binding to nuclear type II sites is 
likely a component of the mechanism of action of this 
bioflavonoid, and type II sites have been identified as 
a component of histone H4, we performed ChIP assays 
with anti-histone H4 antibodies to determine whether 
luteolin modifies the acetylation patterns of H4 associ- 
ated with target gene. We initially evaluated methods for 
enzymatically cleaving MCF-7 cell chromatin (Figure 
6A) to generate approximately 200 bp fragments suit- 
able for the ChIP assays. Subsequently, MCF-7 cells 
were treated with 5 nM estradiol alone (controls) or es- 
tradiol plus 17.5 uM luteolin for 5 hours in phenol-red 
free DMEM supplemented with 10% CSFCS. The cells 
were harvested and chromatin was prepared as in Figure 
6A. Enzymatically digested chromatin from estradiol 
(controls) or estradiol + luteolin treated MCF-7 cells was 
then immunoprecipitated with anti-acetyl histone H4 



Figure 4. GenMAPP Dia- 
gram of Cell Cycle Path- 
way Genes Regulated by 
Luteolin in E2-Treated 
MCF-7 Cells. MCF-7 cells 
were treated with for 5 hours 
with EtOH vehicle, estradiol 
(5 nM), luteolin (17.5 uM) 
or estradiol + luteolin in 2 
uL of EtOH. Three separate 
RNA sets for each group 
were analyzed by Affyme- 
trix HG-U133 PLUS 2.0 
human microarray chips. 
D-Chip output was ana- 
lyzed by GenMAPP to gen- 
erate the diagram. Estradiol 
regulated genes (dotted or 
dashed) were determined by 
comparison of the expres- 
sion in estradiol treated 
groups relative to the EtOH 
controls. Luteolin regulated 
genes were determined by 
comparison of the estradiol 
treated group to the estradiol + 
diol) and significantly (p<0.05) 
results. 
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luteolin treated group. Significantly (p<0.05) upregulated genes are shaded (luteolin) or dotted (estra- 
down-regulated genes are in gray (luteolin) or dashed (estradiol). Table 3 numerically summarizes these 
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antibodies and analyzed by PCR utilizing primers for 
the PLK1 (Figure 6B). Luteolin significantly decreased 
O<0.01) the intensity of the PCR signal for PLK1 (com- 
pare estradiol vs. estradiol + luteolin) in agarose gels. 
Quantitation of the response by scanning the gel bands 
with UN-SCAN-IT software and normalization of the 
signal to the Input DNA's (Figure 6C) confirmed that 
luteolin significantly decreased (/?<0.01) the signal (es- 
tradiol + luteolin) for PLK1 by about 50% relative to 
control (estradiol treated cells). The observations were 
further confirmed by qPCR (Figure 6D) where luteolin 
treatment (estradiol + luteolin) decreased the acetyla- 
tion of histone H4 (K5, K8, K12, K16) associated with 
the PLK1 promoter by 50-60% relative to the estradiol 
controls. These very exciting data suggest that luteolin 
inhibition of PLK1 gene expression (Figures 4 and 5) 
resulted from decreased acetylation of histone H4 asso- 
ciated with the PLK1 gene promoter in MCF-7 cells. We 
believe this response is mediated via luteolin interaction 
with ligand-binding domain on histone H4. 




Figure 5. qPCR Analysis of Cell Cycle Gene Expression. 

MCF-7 Cells were treated as described in Figure 2. RNA from 
these cells was analyzed by real-time PCR (qPCR) as described 
in Figure 3. The estradiol treatment group (used for control) 
represents the 100% value. RNA expression in the estradiol + 
luteolin treatment group is expressed as a % of the estradiol con- 
trol. Data were analyzed statistically by ANOVA and a test on 
the treatment means (Instat, Graphpad Software). */><0.0001. 
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Figure 6. ChIP Assay for 
PLK1 Promoter. Panel 
A. Formaldehyde fixed 
MCF-7 cell chromatin was 
enzymatically digested 
(10-25 minutes) and fol- 
lowing cross-link reversal 
and proteinase K treat- 
ment, the sheared DNA 
was electrophoresed on 1% 
agarose and stained with 
ethidium bromide. The gel 
contains a 100 bp molecu- 
lar weight ladder (MWL), 
non-digested chromatin 
and DNA from 10', 15', 20' 
and 25' digests. The 15'-20' 
digests yielded 100-200 bp 
DNA fragments suitable for 
ChIP assays. Panel B (ChIP 
Assay). DNA cross-linked 
chromatin from estradiol (5 
nM; E, controls) or estradiol 
(5 nM) + luteolin (17.5 uM; 
E + L) treated MCF-7 cells 

was digested to »200 bp (17' digest) and pre-incubated with anti-histone H4 antibodies (Upstate) directed against a peptide consisting 
of amino acid residues 2-19 of histone H4 acetylated at lysines 5, 8, 12 and 16. No signal was observed in the minus Ab controls (not 
shown). Following incubation with protein-A-Sepharose beads, and cross-link reversal, immunoprecipitated DNA was de-proteinized 
and subjected to PCR analysis (33 cycles) with forward and reverse primers for the PLK1 promoter (-262 to -80 bp; (26)). PCR products 
from the input DNA (1:10 dilution) were additional controls. The results from 3 PCR reactions (as in panel B) for estradiol or estradiol + 
luteolin treated MCF-7 cells were quantified with Un-Scan-IT (Silk Scientific Software) and analyzed with Instat (panel C). In Panel D, 
the PLK1 promoter activity in panel B was also quantified by qPCR. Data in the estradiol (controls) or estradiol + luteolin treated samples 
were normalized to the Input DNA for the +Ab samples. Data in panels B-D show that luteolin treatment significantly decreased histone 
H4 acetylation at the PLK1 promoter. 



Chip 


Input 


E E + L 


E E + L 




ChIP E + L Input E 
Sample 



Input E + L 



June 2011 Vol.7 No. 2 Int J Biomed Sci www.ijbs.org 



GENE REGULATION BY LUTEOLIN IN BREAST CANCER CELLS 



DISCUSSION 

The antiestrogenic and antiproliferative activities of 
bioflavonoids such as luteolin in a wide array of biological 
systems are well documented. Compounds like luteolin 
block estrogenic response in classical systems such as the 
rat uterus and also control the proliferation of a wide vari- 
ety of estrogen-independent tissues and cell types includ- 
ing breast, prostate, colorectal, leukemia, meningioma, 
and ovarian cancer cells (1, 2, 7-15). This broad-spectrum 
activity strongly suggests a multiplicity of mechanisms 
are involved in the response profiles of normal and ma- 
lignant cell types to compounds such as luteolin. Stud- 
ies in our laboratory over the past two decades regard- 
ing the interaction of bioflavonoids with nuclear type II 
sites provided a unifying feature regarding their control 
of normal and malignant cell proliferation (1, 2, 8, 16, 18). 
The recent discovery that nuclear type II sites represent a 
ligand-binding domain on histone H4 suggested that type 
II site ligands might regulate gene transcription through 
this binding interaction (3-6, 19). The present microarray 
and qPCR studies have confirmed and extended this con- 
cept and have identified specific genes in the ESP and CCP 
regulated by luteolin which are consistent with the anti- 
estrogenic and antiproliferative properties of this bioflavo- 
noid in numerous experimental systems noted throughout 
this manuscript. It should be noted here that the present 
qPCR studies assessing luteolin effects on ESP or CCP 
genes at 5 hours following treatment were performed only 
to confirm the microarray data identifying these pathways 
as targets for luteolin in MCF-7 cells. Therefore, western 
analyses on protein expression were not conducted at this 
early time since luteolin-induced changes in protein ex- 
pression is typically not observed until approximately 24 
hours following treatment (3, 19). 

One of the most significant observations noted by the 
studies described in this manuscript is the fact that luteolin 
modulated the expression of a number of genes in the ESP, 
and this response certainly is consistent with the antiestro- 
genic properties of luteolin. Even though the compound 
does not bind to either ERa or ER(3, luteolin will block 
estrogenic response (1, 21, 28). Thus, the observation that 
luteolin modulates the expression of a number of genes 
in the ESP supports these earlier studies. Of the luteolin- 
regulated genes in MCF-7 cells (Table 2 and Figures 2 and 
3) all are closely coupled to transcription or cell prolifera- 
tion. TAF9 is a general transcription factor or TATA-box 
binding protein that positions RNA polymerase II on the 
promoter (29). GTF2H2 is a general transcription factor 



(30) and GRB2 is an EGF receptor binding protein (31). 
EGF is mitogenic in breast cancer cells (32). NCOR1 is a 
nuclear co-repressor that mediates ligand-dependent tran- 
scription repression (33) and NCOA3 is a nuclear hormone 
receptor activator that enhances transcriptional activation 
(34). POLR1D is a DNA-directed RNA polymerase (35) 
and RAS is a member of the ras oncogene family involved 
in GTP/GDP-binding and GTPase activity (30). POLR2A 
is an RNA polymerase involved in RNA synthesis (36) 
and DDX5 is an RNA helicase involved in RNA struc- 
ture, splicing and translation initiation (30). Thus, luteo- 
lin regulation of these ESP genes should profoundly af- 
fect estradiol- modulated breast cancer cell proliferation. 
Interestingly, we have recently reported that luteolin also 
blocks EGFR gene transcription and a number of genes in 
the EGFR signaling pathway in PC-3 human prostate can- 
cer cells demonstrating that this bioflavonoid apparently 
targets growth regulatory pathways (3). 

Regular effects of luteolin were also observed on a 
number of genes in the CCP including CCNA2, PLK1, 
PCNA and CDKN1A (Table 3 and Figures 4 and 5). That 
luteolin inhibited CCNA2 expression is also consistent 
with its inhibitory effects on cell proliferation (37). Polo- 
like kinase 1 (PLK1) is involved in entry into and exit 
from mitosis and bipolar spindle formation (38, 39). PCNA 
was upregulated in estradiol stimulated MCF-7 cells and 
down-regulated by combination estradiol plus luteolin 
treated cells (Table 3), suggesting a direct antagonism of 
an estrogen-stimulated gene by luteolin. The remaining 
genes were inhibited by luteolin and not affected by estra- 
diol in MCF-7 cells. PCNA has multiple functions in the 
control of cell proliferation (40) as high levels of PCNA, in 
the absence of p53, are associated with DNA replication. 
High levels of PCNA, in the presence of p53, are associat- 
ed with DNA repair. If PCNA is non-functional or absent, 
apoptosis occurs (40). Thus, luteolin suppression of PCNA 
expression should have a plethora of effects on cell prolif- 
eration depending upon what other p53 related proteins are 
modulated by luteolin. Similarly, luteolin stimulated the 
expression of CDKN1A (Table 3 and Figures 4 and 5), a 
finding consistent with the observation that an increase in 
CDKN1A gene transcription is associated with the inhibi- 
tion of cellular proliferation (41). CDKN1A interacts with 
CDK2 to block cell cycle progression (42). 

On the basis of these observations, we conclude that 
luteolin regulation of breast cancer cell proliferation oc- 
curs through the regulation of genes in both the ESP and 
CCP. Obviously, one has to consider the possibility that the 
effects of luteolin on CCP genes could, in part, result from 
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the down-regulation of ESP genes which control CCP 
genes (43). However we also observed a direct effect of 
luteolin on the expression of the PLK1 gene (Figures 4 and 
5). Luteolin clearly blocked the acetylation of histone H4 
associated with the PLK1 promoter (Figure 6). This was 
also the case for luteolin regulation of the PLK1 gene in 
PC -3 human prostate cancer cells which are not estrogen- 
dependent (19). Taken together, these observations support 
our hypothesis that luteolin regulation of target genes in 
the ESP and CCP may occur through binding interactions 
with type II sites on histone H4 associated with the gene 
promoters. Such binding may affect the recruitment of his- 
tone acetylases or histone deacetylases to these promot- 
ers to modulate the acetylation state of histone lysine tails 
controlling gene transcription (3, 19). We are currently 
evaluating these possibilities in detail. 
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